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ABSTRACT 

We describe the Fortran code CPsuperH2.0, which contains several improvements and 
extensions of its predecessor CPsuperH. It implements improved calculations of the Higgs- 
boson pole masses, notably a full treatment of the 4x4 neutral Higgs propagator ma- 
trix including the Goldstone boson and a more complete treatment of threshold effects in 
self-energies and Yukawa couplings, improved treatments of two-body Higgs decays, some 
important three-body decays, and two-loop Higgs-mediated contributions to electric dipole 
moments. CPsuperH2.0 also implements an integrated treatment of several 5-meson ob- 
servables, including the branching ratios of Bg — > fi'^f^~, Bd — > r"'"r~, Bu — ^ rz/, B — ^ Xg^i 
and the latter's CP- violating asymmetry v4.cp, and the supersymmetric contributions to 
the -B°^ — mass differences. These additions make CPsuperH2.0 an attractive inte- 
grated tool for analyzing supersymmetric CP and flavour physics as well as searches for 
new physics at high-energy colliders such as the Tevatron, LHC and linear colliders. E 



*The program may be obtained from http://www. hep. mEin.ac.uk/u/jslee/CPsuperH. html 
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1 Introduction 



With the imminent advent of the LHC, particle physics experiments are poised to explore 
the TeV energy range directly for the first time. There are several reasons to expect 
new physics in this energy range, such as the origin of particle masses and electroweak 
symmetry breaking, the hierarchy problem and the nature of dark matter. In parallel with 
the direct exploration of the TeV scale, precision experiments at low energies continue to 
place important constraints on the possible flavour and CP-violating structure of any TeV- 
scale physics. Prominent examples include experiments on B and K mesons, and probes of 
electric dipole moments [1]. It is clearly desirable to develop computational tools that can 
be used to calculate consistently observables for both low- and high-energy experiments in 
a coherent numerical framework. This is particularly desirable in view of the possibility 
that the dominance of matter over antimatter in the Universe may be due to CP- violating 
interactions at the TeV scale [2]. 

Supersymmetry is one of the most prominent possibilities for new TeV-scale physics, 
and the minimal supersymmetric extension of the Standard Model (MSSM) provides a 
natural cold dark matter candidate as well as stabilizing the electroweak scale and facil- 
itating the unification of the fundamental interactions. There are many computational 
tools available for calculations within the MSSM. The first to include CP-violating phases 
was CPsuperH [3] based on the renormalization-group-(RG-)improved effective potential 
approach. The Higgs-boson pole-mass shifts are calculated by employing the RG-improved 
diagrammatic approach. The recent versions of FeynHiggs [4] are based the Feynman dia- 
grammatic approach. There are merits in both approaches and the difference between two 
programs may be attributed to some unknown higher-order corrections. 

Some of us have recently published an analysis of several 5-physics observables taking 
into account the most general set of CP-violating parameters allowed under the assumption 
of minimal flavour violation in the supersymmetric sector [5] . For this purpose we used an 
updated and extended computational tool, CPsuperH2 . 0, which we introduce and describe 
in this paper. 

The main new features of CPsuperH2 . are its inclusion of a number of B observables, 
including the branching ratios of Bs Bd t^t~ ^ Bu tu, B Xg'y and 

the latter's CP- violating asymmetry Acp, and the supersymmetric contributions to the 
Bg^ — Bg^ mass differences. In addition, CPsuperH2.0 includes a more complete treatment 
of Higgs-boson pole masses, based on a full treatment of the 4x4 neutral Higgs propagator 
matrix including the Goldstone boson and a more complete treatment of threshold effects in 
self-energies and Yukawa couplings. It also includes improved treatments of two-body Higgs 
decays, some important three-body decays, and two-loop Higgs-mediated contributions to 
electric dipole moments. Therefore, CPsuperH2 . provides an essentially complete, self- 
contained and consistent computational tool for evaluating flavour and CP-violating physics 
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at energies up to the TeV scale. 

The structure of this paper is as follows. Several updated features of CPsuperH2 . are 
described in Section 2. In particular, in Subsection 2.1 we introduce the improved treatment 
of Higgs-boson pole masses, and Section 2.2 contains a description of the improvements in 
the treatment of Higgs decay modes. Then, in Section 3 we describe the CPsuperH2 . 
treatment of two-loop Higgs effects on electric dipole moments. The most important new 
features are described in Section 4, where we discuss its treatment of B observables. In 
each Section, we illustrate in figures some typical results obtained using CPsuperH2 . 0. 

2 Updated Features of CPsuperH2 . 

It is to be understood that, throughout this paper, we follow the notations and conventions 
defined and adopted in CPsuperH for the mixing matrices of neutral Higgs bosons, charginos, 
neutralinos and third-generation sfermions, as well as their masses and couplings, etc. The 
updates to the original version of CPsuperH [3] that are presented here reflect, in part, 
feedback from users, as well as extending it to B observables. 

New common blocks /HC_RAUX/ and /HC_CAUX/ have been introduced for the general 
purpose of storing new numerical outputs which are available in CPsuperH2 . 0: 

• COMMON /HC_RAUX/ RAUX_H 

• COMMON /HC_CAUX/ CAUX_H 

The two arrays RAUX_H and CAUX_H are NAUX=999 dimensional and only parts of them are 
being used presently as shown in Tables [T] and [21 The contents of these two new arrays 
are explained in the corresponding following subsections. These common blocks can also 
be used by users for their specific purposes. 

2.1 Improved Treatment of Higgs-Boson Masses and Propaga- 
tors 

In CPsuperH2 . we make three main improvements in the calculation of the Higgs-boson 
pole masses. 

• The finite threshold corrections induced by the exchanges of gluinos and charginos 
have been included in the top- and bottom-quark self-energies of the neutral and 
charged Higgs bosons. For the explicit expressions of the self-energies, we refer to 
Eqs. (B.14), (B.15), and (B.16) of Ref. [6] 0. 

'''We find that overall minus signs are missing in the expressions of nf^*-22(s)- 
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Figure 1: The masses of the neutral Higgs bosons as functions of tan (3 for the CPX scenario 
[8] taking $3 = = 90° the convention $^ = 0, Msusy = 0.5 TeV, and the charged 

Higgs-boson pole mass Mh± = 160 GeV. In each frame, the dashed line is for the case 
IFLAG_H(12)= 1 and the solid line for other case indicated. 



• Also included are the threshold corrections to the Yukawa couplings \ht^b\ in the one- 
loop running quartic couplings, Xi^\Q = 'rnf"^^) with i = 1 — A. For the explicit 
expressions of \f \ we refer to Eqs. (3.3)-(3.6) of Ref. [7]. 

• An improved iterative method has been employed for the calculation of the pole 
masses. 



As a help in assessing the improvements in the calculation of Higgs sector, new flags 
IFLAG_H(12) and IFLAG_H(60) have been introduced as follows: 



• IFLAG_H(12): 
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Figure 2: The masses of the neutral Higgs bosons as functions of the common phase $a 
for the trimixing scenario [9] taking $3 = —90°. Specifically, in this scenario, tan/3 = 50 
and Mh± = 155 GeV. The lines are the same as in Fig. [H 



— IFLAG_H(12)= 1: Gives the same result as that obtained by the older version of 
CPsuperH. 

~ IFLAG_H(12)= 2: Includes only the threshold corrections to the neutral and 
charged Higgs-boson quark self-energies. 

— IFLAG_H(12)= 3: Includes only the threshold corrections to . 

— IFLAG_H(12)= 4: Includes only the iterative method for the pole masses. 

— IFLAG_H(12)= 5 or 0: All the improvements are fully included. 

• IFLAG_H(60)= 1: This is an error message that appears when the iterative method 
for the pole masses fails. 
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The improvement in the threshold corrections to the top- and bottom-quark Yukawa cou- 
phngs is important when tan/3 is large and the charged Higgs boson is light. In Figs. [1] 
and [2], we show the pole masses of the neutral Higgs bosons for the CPX [8] and trimix- 
ing [9] scenarios, respectively, when IFLAG_H(12) = 2-5 as indicated. In each frame, the old 
calculation with IFLAG_H(12)= 1 (dashed line) is also shown for comparison. 

Finally, RAUX_H(l-6), RAUX_H ( 10-36) , and CAUX_H(l-2) are allocated for numerical 
information on the Higgs-sector calculation based on a renormalization-group-improved 
diagrammatic approach including dominant higher-order logarithmic and threshold correc- 
tions [6,7], see Tables [1] and [2l 
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Figure 3: The absolute value of each component of the neutral Higgs-boson propagator 
matrix D^'^ (s) with (red solid lines) and without (black dashed lines) including off-diagonal 
absorptive parts in the trimixing scenario with $a = —$3 = 90° and IFLAG_H(12)= 5. We 
note that \D^l{s)\ = 1. The three Higgs-boson pole masses are indicated by thin vertical 
lines. 
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In situations where two or more MSSM Higgs bosons contribute simultaneously to a 
process, the transitions between the Higgs-boson mass eigenstates need to be considered 
before their decays. For this reason, we include the complete 4 x 4-dimensional propagator 
matrix (s) spanned by the basis {Hi, H2, H3, G°) [10], including off-diagonal absorptive 
parts [9]. The dimensionless neutral Higgs-boson propagator matrix is given by 

/ s - Mfj^ + i'^mUnis) i'^mUuis) iQmfli3{s) i^mflu{s) \ 

iQmn2i{s) s — Mfj^ + iQmIl22{s) iQmn23{s) iQ'mn24(s) 

iQ^mIl3i(s) i9'mll32(s) s — M]j^ + iQmIl-s3{s) i'^mll34,{s) 

\ i9mn4i(s) i3mll42(s) i3mll43(s) s -|- «Q'mIl44(s) / 



where Mhi 2,3 one-loop Higgs-boson pole masses, and higher-order absorptive effects 

on M/^j 2 3 have been ignored [6]. The label '4' refers to the would-be Goldstone boson of the 
Z boson. The absorptive part of the Higgs-boson propagator matrix receives contributions 
from loops of fermions, vector bosons, associated pairs of Higgs and vector bosons, Higgs- 
boson pairs, and sfermions: 

53m%(s) = 53mnf/ (s) + 53mn5^(s) + Qmflf/{s) + 53mng^(s) + 53mnf/(s) , (2) 

respectively. We refer to Ref. [9] for their explicit expressions. For the Goldstone-Higgs 
mixings, Q'mHi4^4j and SmH44, we take the leading contributions ignoring all gauge-coupling 
mediated parts. We also include the 2 x 2-dimensional propagator matrix for the charged 
Higgs bosons D^^{s) spanned by the basis including off-diagonal absorptive 

parts: 

The relevant Goldstone-boson couplings are given in Appendix [Bl For the 16 elements of 
the neutral Higgs-boson propagator matrix D^°{s) and for the 4 elements of the charged 
Higgs-boson propagator matrix D^'^{s), the slots CAUX_H(100-119) are used as shown in 
Table [21 In Fig. [3l as an example, we show the absolute value of all components of the 
Higgs-boson propagator matrix (s) as functions of for the trimixing scenario with 
$^ = -$3 = 90°. 

It is important to remark that the 4x4 propagator matrix ([1]) is sufficient to encode all 
Hi—Z- and G^—Z mixing effects within the Pinch Technique (PT) framework [10,11], which 
has been adopted here to remove consistently gauge-dependent and high-energy unitarity- 
violating terms from 53mHjj(s) [9]. For example, the self-energy transition Hi Z^, 
^ZH, = P^^ZHi, is related to Hgoh, through 

sUzhAS) = -tMlUcOHAS) , (4) 



with s = p^. We recall that the self-energy transitions Hi and G*^ ^ 7 are completely 
absent within the PT framework. More details may be found in [10]. 

Note that the elements of the propagator matrix depend on the center-of-mass energy, 
denoted by a/S, which is stored in RAUX_H(101), see Tabled] Along with D^°'^'^{s), the 
s-dependent couplings of the neutral Higgs bosons to two gluons, S'f(V^) and Pf(V^), 
and two photons, Sj{\/§) and P^{^/I), are needed when we consider the production of 
the neutral Higgs bosons and study its CP properties at the LHC [9, 12, 13] and a 77 
colhder [14-16]. They are calculated and stored in CAUX_H(130-135) and CAUX_H(140-145) 
as shown in Table [2l We have included the dominant contributions coming from the 
tan f3 enhanced loops of sbottoms and gluinos and the subdominant ones coming from the 
stop-higgsino mediated diagrams. Also included are the resummed corrections to Yukawa 
couplings. For the electroweak corrections, see next subsection. For the next-to-leading- 
order QCD corrections, appropriately calculated K factors should be taken into account 
separately in the calculation of production cross sections [17, 18]. 

Two additional flags are used to control the inclusion of the off-diagonal absorptive 
parts and print out the the s-dependent propagator matrix and the s-dependent Higgs 
couplings to two photons and gluons: 

• IFLAG_H(13)= 1: Does not include the off-diagonal absorptive parts in the propagator 
matrices '^^{s). 

• IFLAG_H(14)= 1: Prints out each component of the Higgs-boson propagator matrices 
D^'^'^^{s) and the s-dependent couplings ^/'^(V^) and Pi''^{VS). 



2.2 Improved Treatment of Higgs-Boson Couplings and Decays 

The main updates include: 

• The electroweak corrections to the neutral Higgs couplings to pairs of tau leptons and 
6-quarks [19]. The explicit formulae used in the code for the corrections, including 
non- vanishing CP phase effects, could be found in Ref. [9] and Eqs.(A.l)-(A.2) of 
Ref. [3]. 

• The three-body decay t*b — > W~^bb. Some three-body decays play important 
role in Higgs searches [20]. In addition to the three-body decays involving more than 
one massive gauge boson considered previously, we include the three-body decay 

— > t*b W'^bb in the new version. The decay width is given by 

T{H+ W^bb) = 

5127r'^ Jo [i- - X2 - Kt + Kb) + l^tlt 
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where = ■ml/Mfj±, 74 = Tf/Mfj± and Xi = 2Ei/MH± with Ei and E2 being the 
energies of the b and b quarks, respectively. In the charged Higgs-boson rest frame, 
the function E{xi,X2) is given by 



F{xi,X2) 



\9l\ 



+ \9r\ 



2xi + 2x2 — 3 + 2kw — 2Ki,i%t 



" (1 -X2) 
xl + X1X2 — Sxl — 2xiX2 + 3X2 + Xi — 1 
+ {x\ + 2xiX2 — 4x2 — 2xi + 3 — 2kw) 



+2^jKbKt ^e{gLg*R 



-Kb ^ — 2xi + 3 + 2kw 

(X2 - P2 



-2x2 ~ X1X2 + 5x2 + xi — 3' 

Kw 



Kw 



+ {—X2 + 1 — 2kw) + 2Kb 



2Kt 



(6) 



where ql = gfj+jb - igu+tb and Qr = gfj+jf^ + ig^+^f,. 

The contributions from tau-lepton and charm-quark loops to the couplings S2{Mh.) 
and P^{Mh:). 

A new flag IFLAG_H(57)= 1: This is an error message that appears when one of the 
magnitudes of the complex input parameters is negative. 



The CPsuperH homepage has been continuously brought up to date after its first appearance 
to include the updates discussed in this subsection and others not mentioned here. We refer 
to the file 0LIST_V1 for a full list of updates to the original version which can be found in 
the CPsuperH homepage. 



3 Higgs-Mediated Two-Loop Electric Dipole Moments 

The CP phases in the MSSM are significantly constrained by measurements of Electric 
Dipole Moments (EDMs). In particular, the EDM of the Thallium atom may provide 
currently the most stringent constraint on MSSM scenarios with explicit CP violation. 
The atomic EDM of ^os^jj gg^g j^g ^lain contributions from two terms [21,22]: 

dT\ [e cm] = -585 ■ 4 [e cm] - 8.5 x 10"^^ [e cm] ■ {Cs TeV^) + ■ ■ ■ , 

= {dTiT [e cm] + (^Ti)^^ [e cm] + ■ ■ ■ , (7) 

where de denotes the electron EDM and Cs is the coefficient of the CP-odd electron- 
nucleon interaction Ccs = ei'-f^e NN. The dots denote sub-dominant contributions 
from 6-dimensional tensor and higher-dimensional operators. 
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Figure 4: The Thallium EDM d-n = d^i x 10^^ [ecm] in the CPX scenario with $^ = 
$3 = 90° and Msusy = 0.5 TeV taking IFLAG_H(12)= 5 [25]. The different shaded regions 
correspond to different ranges of |(iTi| o,s shown. Specially, in the narrow region denoted by 
black squares, one has < 1; consistent with the current Thallium EDM constraint. 



The contributions of the first- and second-generation phases, ^ ^^'i s ' EDMs 
can be drastically reduced either by assuming that these phases sufficiently small, or if the 
first- and second-generation squarks and sleptons are sufficiently heavy. However, even 
when the contributions of the first and second generation phases to EDMs are suppressed, 
there are still sizeable contributions to EDMs from Higgs-mediated two- loop diagrams [23]. 

The Higgs-mediated two-loop Thallium (^^1)5 electron (d^), and muon (d^) EDMs 
are calculated and stored in RAUX_H(111-120) as shown in Table [H The Thallium and 
electron EDMs consist of: 

= idfY + {d^f+idfy + {dff + {dfr\ (8) 
The explicit expressions for the EDMs in the CPsuperH conventions and notations may be 
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found in Ref. [24]. A flag IFLAG_H(15)= 1 is used to print out the results of the EDM 
calculations: 

• IFLAG_H(15)= 1: Print out EDMs. 

In Fig. m we show the rescaled Thallium EDM din = x 10^^ in units of ecm 
in the tan jS-Mn^ plane, in the CPX scenario with IFLAG_H(12)= 5. We observe, when 
tan/3 < 5 and Mh-^ ^10 GeV, one may have \dTi\ < 1 only in the narrow region denoted 
by black squares which is consistent with the current 2-a upper bound on the Thallium 
EDM [26]: [dxil < 1-3 x lO^^'^ [ecm]. We note that the region 8 GeV < Mh, < 10 
GeV with tan/3 < 10 has not been excluded by the combined constraints from the LEP 
searches [27] and the T(15') jHi decay [28]. 

The Thallium EDM constraint can be evaded by assuming cancellations between the 
two-loop contributions considered here and possible one-loop contributions which depend 
on different CP-odd phases related to the first and second generations of squarks and 
sleptons. For example, assuming cancellation of less than 1 part in 10, the region with 
1 < |(ixi| < 10 in Fig. H] is allowed. In the future, this treatment of the most important 
two-loop contributions to the Thallium EDM will be supplemented by a more complete 
implementation of calculations of the well-known 1-loop contributions to this and other 
EDMs. 

4 5-Meson Observables 

An important innovation in CPsuperH2 . is the inclusion of the following important Higgs- 
mediated S-meson observables: 

• The branching ratio of Bg meson into a pair of muons: B{Bs fifi), 

• The branching ratio of Bd meson into a pair of tau leptons: B{Bd tt), 

• The SUSY contribution to the B^-B^ mass difference: AM^^^ , 

• The SUSY contribution to the B^-B^ mass difference: AM^^^ , 

• The ratio of the branching ratio B{Bu tv) to the SM value: 



• The branching ratio B{B ^ ^sl) and the direct CP asymmetry Acp{B ^s7). 



Rb 



•tu — 



B{B~ 
B^^iB, 



u 
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We adopt the most recent gauge-invariant and flavour-covariant formahsm to calculate 
the flavour- changing effective Lagrangian for the interactions of the neutral and charged 
Higgs fields to the up- and down-type quarks including a new class of dominant subleading 
contributions [5]. In the current version, the single-Higgs insertion approximation is used. 

For the calculations of 5-meson observables, the array SMPARA_H for the SM param- 
eters has been extended to include information on the CKM matrix, parameterized via A, 
A, p, and f/, as seen in Table [31 The CKM matrix is constructed as [29] 



V 



Cl2 Ci3 



— Sl2 C23 — C12 S23 Si3 e 
V S12S23 - Ci2C23Si3e*'^ 



iS 



S12 Cl3 
C12 C23 — ■5i2 S23 Si3 e 



iS 



-C12S23 - Si2C23Si3e' 



S23 Cl3 
C23 Cl3 / 



(9) 



where Sij = sin^jj, Cij = cosOij, and 6 is the KM phase with Sij , Cij > 0. In terms of A, A, 



p, and T], they are given by 

A , Sos = AX^ 



S12 



si3e' 



A\^{p + ifi)Vl- A^X^ 



(10) 



and Cij = yl — |sijP. The SUSY parameter array SSPARA_H is also extended to include the 
hierarchy factors Pq^^dXe between the first two and third generations [30], see Table HI 
In the super-CKM basis, the 3x3 squark mass matrices squared are taken to be diagonal: 





= "^la 


X 


diag (p|,p|,l). 




= "^k 


X 


diag(p|,p|, 1) , 


ml 


= "^X 


X 


diag {p%pI,1), 


Mi 


= 


X 


diag (p|, pI, 1) , 


Ml 


= 


X 


diag (p|,p|,l). 



:iii 



Finally, the results for the i?-meson observables are stored in RAUX_H(130-136) as 
shown in Table [TJ The SUSY contributions to the AB = 2 transition amplitudes are 
stored in CAUX_H(150) and CAUX_H(151), see Tabled Note the relations RAUX_H(132) = 
2 X |CAUX_H(150)| and RAUX_H(133) = 2 x |CAUX_H(151) |. Two fiags IFLAG_H(16) and 
IFLAG_H(17) are used to print out the results of the calculation of 5-meson observables: 



• IFLAG_H(16)= 1: Print out 5-meson observables. 

• IFLAG_H(17)= 1: Print out details of the B Xs'y calculation. 



For numerical examples of 5-meson observables, we take the CPX scenario [8] with 
^susY = 0.5 TeV and the common A-term phase $a = ^At = ^At = ^a^ iii the convention 



12 



$^ = 0°. We take account of the dependence on the hierarchy factors Pqfi^o between the 
first two and the third generations, taking a common value p for the three of them. 

Figure [5] shows the dependence of the branching ratio B{Bs — * p^ p~) on the phase 
of the gluino mass parameter $3 for four values of tan/3. The charged Higgs-boson pole 
mass is fixed at Mh± = 200 GeV. In each frame, two sets of three lines are shown. The 
upper lines are for higher p = 10 and the lower ones for p = 1. For fixed p, three lines show 
the cases of $^ = 0° (solid), 90° (dashed), and 180° (dash-dotted). The p dependence is 
shown in Fig. [6l We clearly see the GIM operative point mechanism discussed in Ref. [30] 
around p ~ 1.2 when (<I>3 , $a) = (0° , 180°) (solid lines). Figure [7] shows the rescaled 
branching ratio = B{Bs P^P~) x 10'' in the Mj^j-tan/? plane when the phases are 
fixed at $^ = $3 = 90°. The unshaded region is not theoretically allowed. Only the 
region with B^ < 0.58 is consistent with the current experimental upper limit at 95 % 
C.L., corresponding to tan/5 < 20 (8) for p = 1 (10). 

The rescaled branching ratio Bsj = B{B ^ ^sl) x 10^ is shown in Fig. [HI In 
contrast to the Bg — * p^ p~ case, we observe that higher tan/? region is experimentally 
allowed: tan/5 > 35(20) for p = 1 (10). This is because the charged-Higgs contribution 
is suppressed due to the threshold corrections when tan/5 is large. The charged-Higgs 
contribution io B ^ is proportional to 1/(1 + |fi:ptan^/5) [33], where n represents the 
threshold corrections with |k| ~ 0.05 for the parameters chosen [34]. 

Figure [9] shows the ratio of the branching ratio B{Bu tu) to its SM value, Rbtv 
In the left frame with p = 1, we see two connected bands of the experimentally allowed 
1-(T region, 0.62 < Rbtv < 1-38. If we consider the 2-a limit, only the upper-left region 
with ^ 95 GeV and tan/5 > 35 is not allowed. For larger p = 10, the allowed region 
becomes narrower. 

In Fig. [ini we show the region satisfying the experimental constraints from B{Bs 
p^p~) (95 %), B{B — > Xs7) (2 a), and Rbtu (1 cr)- First we observe that there is no 
region that satisfies the Bg P~^P~ and B —>■ Xs'~f constraints simultaneously for both 
p = 1 and 10. If one neglects the constraint from B{Bs p~^p~), only the high-tan /5 
region would remain. Taking account of B^ ru constraint, the region with tan/5 > 36 
and M//^ > 80 GeV is allowed when p = 1. On the other hand, neglecting the constraint 
from B{B ^ ^sl), the allowed region is constrained in the parameter space with tan /5 < 20 
and > 10 GeV for p = 1. For p = 10, the B — > Xg'y constraint is relaxed but those 
from B{Bs P-'^P') and Rbtu become more stringent. 

Finally, in Fig. [TTl we show the region allowed experimentally by the measurement 
B{Bd T^T~) < 4.1 X 10~^ (90 %) [36] (upper frames) and the regions where the SUSY 
contribution is smaller than the measured values of B^-B^ mass difference [37] (middle 
frames) and B^-B^ mass difference [29] (lower frames). We see that the B{Bct t^t~) 
constraint has the least impact on these parameter planes, whereas the impacts of the 
Bg-Bg and B^-B^ mass differences are similar. 
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These examples illustrate the possible interplays between the different 5-meson ob- 
servables, and how they may vary significantly with the values of the CP- violating phases. 
CPsuperH2.0 provides a unique tool for combining these constraints and pursuing their 
implications for other observables. In the future, the CPsuperH2 . treatment of these im- 
portant i?-meson observables will be supplemented by the implementation of calculations 
of other flavour observables, including the K sector. 

5 Summary and Outlook 

We have presented in this paper a description of the new features of the Fortran code 
CPsuperH2 . 0. In addition to improved calculations of the Higgs-boson poles masses with 
more complete treatment of threshold effects in self-energies and Yukawa couplings, the 
complete A^A (2x2) neutral (charged) Higgs-boson propagator matrices with the Goldstone- 
Higgs mixing effects have been consistently implemented. Specifically, the neutral Higgs- 
boson propagator matrix constitutes a necessary ingredient for the studies of a system 
of strongly-mixed Higgs bosons at colliders together with the center-of-mass dependent 
Higgs-boson couplings to gluons and photons. It also provides the improved Higgs-boson 
couplings to tau leptons, b quarks, and two photons. The important three-body decay 
H+ t*h W+bb is included. 

In order to provide a more complete, consistent tool for calculating CP-violating 
observables in the MSSM, and specifically to incorporate the important constraints coming 
from precision experiments at low energies, CPsuperH2.0 has been extended to include 
a number of i?-meson observables, as well as the Higgs-mediated two-loop contributions 
to EDMs of the Thallium atom, electron and muon. The currently available 5-meson 
observables are the branching ratios of Bg yU'^/U^, Bd t~^t~ , Bu — > tu, B — > and 
the latter's CP- violating asymmetry ^cP) and the supersymmetric contributions to the 
B^^ — B^^ mass differences. Further low-energy observables are to be included in future 
updates. 

The improved Fortran code CPsuperH2 . provides a coherent and complete numerical 
framework in which one can calculate consistently observables in both low- and high-energy 
experiments probing physics beyond the SM. 
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A List of changes 



Here we summarize the improved features introduced in CPsuperH2 . compared to the 
prior version of CPsuperH. 

• New common blocks: 

— COMMON /HC_RAUX/ RAUX_H(NAUX=999) , see Table [1] 

— COMMON /HC_CAUX/ CAUX_H(NAUX=999) , see Table [2] 

• Extended arrays for input parameters: 

— SMPARA_H(NSMIN=19), see Table [3] 

— SMPARA_H(NSSIN=26), see Table H 

• New names for improved FORTRAN files: 

— cpsuperh.f — > cpsuperh2.f 

— fillpara.f — > fillpara2.f 

— fillhiggs.f — > fillhiggs2.f 

— fillcoupl.f — > fillcoupl2.f 

— fillgambr.f — > f illgambr2 . f 

• New FORTRAN files: 

— filldhpg.f is to calculate the full propagator matrices D^°'^^(s) and the s- 
dependent couplings Sf''^{\^) and Pf ''^ {s/I) . 

— higgsedm.f is to calculate Higgs-mediated two-loop EDMs of Thallium, elec- 
tron, and muon. 

— f illbobs . f is to calculate the S-meson observables: B{Bs fifi), B{Bd tt), 
AMiysY, AMiysY^ ^^^^^ ^ X,7), and Acp{B ^ X,^). 

• New flags: 

— IFLAG_H(12) =0 — 5: For the level of improvement in the calculation of the 
Higgs-boson pole masses. 

— IFLAG_H(13) = 1: Not to include the off-diagonal absorptive parts in the propa- 
gator matrices '^^{s). 



IFLAG_H(14) = 1: Print out the the elements of the full propagator matrices 
D^'^''^^(s) and the s-dependent couplings Sf'^{VI) and Pf''^{V§). 
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— IFLAG_H(15) = 1: Print out EDMs. 

— IFLAG_H(16) = 1: Print out 5- meson observables. 

— IFLAG_H(17) = 1: Print out B ^ Xs'j details. 

— IFLAG_H(57) = 1: This is an error message that appears when one of the mag- 
nitudes of the complex SUSY input parameters is negative. 

— IFLAG_H(60) = 1: This is an error message that appears when the iterative 
method for the neutral Higgs-boson pole masses fails. 



B Goldstone-boson couplings to third-generation fermions 
and sfermions 

Here we present the Goldstone-(s)fermion-(s)fermion couplings in the CPsuperH convention. 

• G'-f-f 

f=t,b,T ^^^^W 

where 

dcm = -1 , 9G%b = 9G0fr = +1 • (B.2) 

• G^-f-f 

^G^ff = ^ G^ f, {nif^ P, - mf^Pn)h + h.c. (B.3) 

= -gtbG^t{gl;+ji, + ig^+ti,-f5)b- g^^rG-^iyr{gG+p^r + ^9G+u^Tl5)T + h.c. , 



where 



_ grrit g _ 1 - mb/rrit p _ . 1 + rrib/mt 

_ g^T 5_i p _ - ^ (Ti A\ 



G'-f*-f 

f = t,b,T 



G0/7 = - E ^7go/*/,(G°/;/,), (B.5) 



where 
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The couplings in the weak-interaction basis are given by 

V2\ -thtisfsAt-cpfi*) 

_]_( -ihlicpAl-s^fx) 

v/2 V ^ hicpAb - Sf3fL*) 

pGOff ^ J_f -ih*{cf3Al-Sf3fi) 

V2 V ihricfsAr - Sp^*) 



(B.7) 



where 



/:g±/7' = V gaHfb, (G^ ^1 h) + ^' (C^ u* n) + h.c. , (B.8) 

V g^^^^ = (r^^**^) f/i:f/^^. and ^ = T^^^^- t/^, . (B.9) 

The couphngs in the weak-interaction basis are given by 

\ K {spAu - Cfsjj.*) J 

YG^Kr^ = -/i; (c^a; - s^/x) . (B.io) 



C Sample new outputs 

Here we show the new outputs of CPsuperH2.0 for the CPX scenario with tan/3 = 5, 
Mh± = 300 GeV, MgusY = 500 GeV, and $^ = $3 = 90°. 

• IFLAG_H(l) = 1: In the new version, we are using mb(mf°''^) = 3.155 GeV and 
mc{m^°^^) = 0.735 GeV as defaults. Note also that the list of the SM and SUSY 
input parameters is extended to include the CKM matrix and the diagonal sfermion 
mass matrices. 



Standard Model Parameters in /HC_SMPARA/ 



AEM_H = 0.7812E-02 

ASMZ_H = 0.1185E+00 

MZ_H = 0.9119E+02 

SW_H = 0.4808E+00 

ME_H = 0.5000E-03 



alpha_em(MZ) 

alpha_s (MZ) 

Z boson mass in GeV 

sinTheta_W 

electron mass in GeV 
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MMU_H 







1065E+00 


muon mass in GeV 






MTAU_H 








1777E+01 


tau mass in GeV 






MDMT_H 







4000E-02 


d- quark mass at M_t"pole 


in 


GeV 


MSMT_H 







9000E-01 


s-quark mass at M_t"pole 


in 


GeV 


MBMT_H 








3155E+01 


b-quark mass at M_t"pole 


in 


GeV 


MUMT_H 








2000E-02 


u-quark mass at M_t"pole 


in 


GeV 


MCMT_H 








7350E+00 


c-quark mass at M_t"pole 


in 


GeV 


MTPOLE_H 







1743E+03 


t-quark pole mass in GeV 






GAMW_H 








2118E+01 


Gam_W in GeV 






GAMZ_H 








2495E+01 


Gain_Z in GeV 






EEM_H 


_ 





3133E+00 


e = (4*pi*alpha_em)"l/2 






ASMT_H 








1084E+00 


alpha_s (M_t"pole) 






CW_H 








8768E+00 


cosTheta_W 






TW_H 







5483E+00 


taiiTheta_W 






MW_H 







7996E+02 


W boson mass MW = MZ*CW 






GW_H 







6517E+00 


SU(2) gauge coupling gw= 


=e/ 


s_W 


GP_H 







3573E+00 


U(1)_Y gauge coupling gp= 


=e/c_W 


V_H 







2454E+03 


V = 2 MW / gw 






GF_H 







1174E-04 


GF=sqrt(2)*gw~2/8 MW"2 in GeV"-2 


MTMT_H 







1666E+03 


t-quark mass at M_t"pole 


in 


GeV 



CKM Matrix 



1 V_ud 1 




1(0 


9738E+00 





OOOOE+00) 1 







9738E+00 


1 V_us 1 




1(0 


2272E+00 





OOOOE+00) 1 







2272E+00 


1 V_ub 1 




1(0 


2174E-02 




3349E-02) 1 







3993E-02 


1 V_cd 1 




l(- 


2271E+00 




1377E-03) 1 







2271E+00 


|V_cs| 




1(0 


9730E+00 




3213E-04) 1 







9730E+00 


1 V_cb 1 




1(0 


4222E-01 





OOOOE+00) 1 







4222E-01 


1 V_td 1 




1(0 


7478E-02 




3259E-02) 1 







8157E-02 


|V_ts| 




l(- 


4161E-01 




7602E-03) 1 







4162E-01 


1 V_tb 1 




1(0 


9991E+00 





OOOOE+00) 1 







9991E+00 



Real SUSY Parameters in /HC_RSUSYPARA/ 



TB_H 


= 


5000E+01 


tan (beta) 






CB_H 


= 


1961E+00 


cos (beta) 






SB_H 


= 


9806E+00 


sin(beta) 






MQ3_H 


= 


5000E+03 


M_tildeQ_3 


in 


GeV 


MU3_H 


= 


5000E+03 


M_tildeU_3 


in 


GeV 


MD3_H 


= 


5000E+03 


M_tildeD_3 


in 


GeV 
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ML3_H = 0.5000E+03 : M_tildeL_3 in GeV 
ME3_H = 0.5000E+03 : M_tildeE_3 in GeV 



Complex SUSY Parameters in /HC_CSUSYPARA/ 



1 MU_H 1 




0. 


.2000E+04: 


:Mag. 


of 


MU 


parameter 


in 


GeV 


|M1_H| 


- 


0. 


.5000E+02: 


:Mag. 


of 


Ml 


parameter 


in 


GeV 


1 M2_H 1 


= 


0. 


.lOOOE+03: 


:Mag. 


of 


M2 


parameter 


in 


GeV 


1 M3_H 1 




0. 


. lOOOE+04: 


:Mag. 


of 


M3 


parameter 


in 


GeV 


1 AT_H 1 




0. 


. lOOOE+04: 


:Mag. 


of 


AT 


parameter 


in 


GeV 


1 AB_H 1 




0. 


. lOOOE+04: 


:Mag. 


of 


AB 


parameter 


in 


GeV 


1 ATAU_H 1 




0. 


. lOOOE+04: 


:Mag. 


of 


ATAU parameter : 


in GeV 


ARG(MU_H) 




0. 


.OOOOE+00: 


:Arg. 


of 


MU 


parameter 


in 


Degree 


ARG(M1_H) 




0. 


.OOOOE+00: 


:Arg. 


of 


Ml 


parameter 


in 


Degree 


ARG(M2_H) 




0. 


.OOOOE+00: 


:Arg. 


of 


M2 


parameter 


in 


Degree 


ARG (M3_H) 




0. 


.9000E+02: 


:Arg. 


of 


M3 


parameter 


in 


Degree 


ARG(AT_H) 




0. 


.9000E+02: 


:Arg. 


of 


AT 


parameter 


in 


Degree 


ARG(AB_H) 




0. 


.9000E+02: 


:Arg. 


of 


AB 


parameter 


in 


Degree 


ARG(ATAU_H) 




0. 


.9000E+02: 


:Arg. 


of 


ATAU parameter : 


In Degree 



Diagonal Sfermion Mass Matrices [GeV] (Not squared) : 
M_q = 0.5000E+03 X Diag(0 . lOOOE+01 O.lOOOE+01 O.lOOOE+01) 
M_U = 0.5000E+03 x Diag(0 . lOOOE+01 O.lOOOE+01 O.IOOOE+Ol) 
M_D = 0.5000E+03 x Diag(0 . lOOOE+01 O.lOOOE+01 O.IOOOE+Ol) 
M_L = 0.5000E+03 x Diag(0 . lOOOE+01 O.lOOOE+01 O.IOOOE+Ol) 
M_E = 0.5000E+03 x Diag(0 . lOOOE+01 O.lOOOE+01 O.IOOOE+Ol) 



Charged Higgs boson pole mass : 0.3000E+03 GeV 



• 1FLAG_H(2) = 1: The masses and mixing matrix of the neutral Higgs boson change 
due to the improvement in their calculations and the new input for the 6-quark mass. 



Masses and Mixing Matrix of Higgs bosons : 

HMASS_H(I) and OMIX_H(A,I) 



HI Pole Mass = 0.1193E+03 GeV 

H2 Pole Mass = 0.2718E+03 GeV 

H3 Pole Mass = 0.2983E+03 GeV 
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Charged Higgs Pole Mass = 0.3000E+03 GeV [SSPARA_H(2)] 

[HI] [H2] [H3] 

[phi_l] / 0.2457E+00 0.3360E+00 0.9093E+00 \ 
0(IA,IH)= [phi_2] I 0.9693E+00 -.7551E-01 -.2340E+00 I 

[ a ] \ -.9973E-02 0.9388E+00 -.3442E+00 / 



• IFLAG_H(14) = 1: The elements of the propagator matrices (s) and the s- 

dependent couphngs of the neutral Higgs bosons to two photons, S]{\fi) and P^lV^); 
and two gluons, Sf{\f^) and Pf{\/3), taking = Mh2- The couplings are compared 
to their values at the Higgs-boson pole masses: S'7(V^ = Mm) = NHC_H(88, IH), 
P^{^ = Mih) = NHC_H(89,IH), Sf{v^ = M^) = NHC_H(84, IH), Pf{VI = Mjh) = 
NHC_H(85, IH). 



DNH4 at sqrts = 0.2718E+03 GeV 



DNH4 [HI, HI] 




(0 


1238E+01 





2290E- 


-01) 







1238E-H01 


DNH4 [H2,H2] 




(0 


5542E-01 




1611E+04) 







1611E-H04 


DNH4 [H3,H3] 




(- 


4876E+01 




2128E- 


-01) 







4876E-H01 


DNH4 [H1,H2] 




(- 


1607E+00 




2973E- 


-02) 







1608E+00 


DNH4 [H1,H3] 




(- 


5956E-05 





2606E- 


-03) 







2606E-03 


DNH4 [H2,H3] 




(- 


4893E-01 




2377E- 


-03) 







4893E-01 


DNH4 [GO, HI] 




(0 


3403E-06 




1825E- 


-04) 







1825E-04 


DNH4 [G0,H2] 




(0 


1872E+00 





2446E- 


-05) 







1872E-H00 


DNH4 [G0,H3] 




(- 


2222E-06 





5181E- 


-04) 







5182E-04 


DNH4 [GO, GO] 




(0 


lOOOE+01 




2365E- 


-05) 







lOOOE-HOl 



DCH2 at sqrts = 0.2718E-I-03 GeV 



DCH2 [H-H , H-h] : I ( - . 4576E-I-0 1 

DCH2[H+,G-h] : |(-.1256E-03 

DCH2[G+,H-h] : |(-.1256E-03 

DCH2[G-H,G-^] : I (0 . lOOOE-i-Ol 



-.2790E-01)| = 0.4576E-H01 

0.2294E-01)| = 0.2294E-01 

0.2294E-01)| = 0.2294E-01 

-.6202E-03)| = O.lOOOE-HOl 



Comparisons of the H-photon-photon couplings at MH"pole 
and those at sqrt{s} = 0.2718E-I-03 GeV 



S couplings P couplings 

HIPP(M): (-.6615E-I-01 0.6386E-01) (0.1303E-01 0.7314E-03) 
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HIPP(S): (-.3180E+01 -.6078E+01) (0.1779E-01 0.2017E-02) 

H2PP(M): (-.9852E+00 0.3333E-01) (-.6867E+00 -.2221E+00) 

H2PP(S): (-.9852E+00 0.3333E-01) (-.6867E+00 -.2221E+00) 

H3PP(M): (-.4272E+00 0.2509E+00) (0.5178E+00 0.7028E-01) 

H3PP(S): (-.3695E+00 0.2852E+00) (0.4567E+00 0.7475E-01) 



Comparisons of the H-glue-glue couplings at MH"pole 
and those at sqrt{s} = 0.2718E+03 GeV 



S couplings P couplings 



HIGG(M) 


(0 


5792E+00 





4164E-01) 


(0 


5316E-02 




6809E- 


-03) 


HIGG(S) 


(0 


7358E+00 





8932E-02) 


(0 


6510E-02 




1457E- 


-03) 


H2GG(M) 


(- 


3557E+00 





2591E-02) 


(- 


1970E+00 




3456E- 


-01) 


H2GG(S) 


(- 


3557E+00 





2591E-02) 


(- 


1970E+00 




3456E- 


-01) 


H3GG(M) 


(- 


2240E+00 





2860E-01) 


(0 


1855E+00 





2231E- 


-02) 


H3GG(S) 


(- 


2150E+00 





3413E-01) 


(0 


1585E+00 





2662E- 


-02) 



• IFLAG_H(15) = 1: The Higgs- mediated two-loop Thallium, electron, and muon EDMs. 
For the Thallium case, the two main contributions from the electron EDM and the 
CP-odd electron-nucleon interaction are shown separately. 



Higgs-mediated two-loop EDMs 
Phi_3 = 0.9000E-H02"o and Phi_At = . 9000E-^02"o 



Thallium [10"-24 ecm] : -.2612E-I-01 

[-.2568E-I-01 from electron EDM] 
[-.4467E-01 from C_S EDM] 

Electron [10"-26 ecm]: 0.4389E-I-00 

Muon[10"-24 ecm] : 0.8997E-I-00 



• IFLAG_H(16) = 1: The 5-meson observables. 



B Observables 



B(B_s -> mu 
B(B -> X_s 
B(B_u -> tau 



mu ) X 10"7 
gamma) x 10"4 
nu)/B(SM) 



0.3710E-01 
0.4396E-1-01 
. 9854E-I-00 
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B(B_d -> tau tau) x 10"7 = 0.2294E+00 

ACP(B -> X_s gamma) x 10"2 = -.7954E-01 [%] 

Delta M [B_d] (SUSY) = 0.6659E-04 [1/ps] 

Delta M [B_s] (SUSY) = 0.1982E-01 [1/ps] 



• IFLAG_H(17) = 1: The details of the B Xg'-f calculation. As a default, we use 
'nT'cifJ'c = rn^°^^) to capture a part of NNLO corrections [38]. The case when only the 
charged-Higgs contribution is added to the SM prediction is also shown. 



B -> X_s gamma 
delta and E_gamma" cut [GeV] : 0.3333E+00 0.1601E+01 



b-q masses [GeV] (pole, (§mb"pole, (§mt"pole) : 

0.4802E+01 0.4415E+01 0.3155E+01 

c-q masses [GeV] (pole, (9mc"pole, (§mb~pole) : 

0.1415E+01 0.1250E+01 0.1029E+01 

mu_b and mu_c [GeV] : 0.4802E+01 0.1415+01 



BR X 10"4: 0.4396E+01 (SM+Charged Higgs+Chargino) 

[0.4471E+01 (SM+Charged Higgs)] 

[0.3351E+01 (SM)] 

AGP X 10"2: -.7954E-01 % 
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Table 1: The contents of the array RAUX_H. In RAUX_H(22) and RAUX_H(23), the notation 
h^j is for the Yukawa couplings without including the threshold corrections. The notations 
which are not explained in the text follow the conventions o/CPsuperH [3] and Refs. [5-7]. 
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pole 
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Table 2: The contents of the array CAUX_H. The notations which are not explained in the 
text follow the CPsuperH [3] convention. 
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Table 3: The contents of the extended SMPARA_H(IP). 
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Figure 5: The branching ratio B{Bs — > fJ''^fJ'~) x 10'^ as a function of $3 for four values 
of tan /3: tan/5 = 10 (upper left), 20 (upper right), 30 (lower left), and 4O (lower right). 
The CPX scenario is taken with Msusy = 0.5 TeV and Mh± = 200 GeV in the convention 
$^ = 0. In each frame, the lower three lines are for the case p = pg = p^y = pfj = 1 and 
the upper lines for p = 10 where the solid, dashed, and dash-dotted lines are for $^ = 0°, 
90°, and 180°, respectively. The current 95 % experimental upper hound, 5.8 x 10~® [31], 
is also shown as a horizontal line in each frame. 



30 




Figure 6: The branching ratio B{Bs — s> x 10^ as a function of the common hierarchy 

factor p ^ Pq = pjj = pj^ for four values of tan P: tan/5 = 10 (upper left), 20 (upper 
right), 30, (lower left), and 40 (lower right). The CPX scenario is taken with Msusy = 0.5 
TeV and Mh± = 200 GeV in the convention $^ = 0. In each frame, the solid line is for 
($3,$^) = (0° , 180°) and the dashed one for (90°, 90°). The current 95 % experimental 
upper hound, 5.8 x 10~® [31], is also shown as a horizontal line in each frame. 
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Figure 7: The branching ratio = 
The CPX scenario is taken with $^ 
of the common hierarchy factor: p = 
theoretically allowed. The different shaded regions correspond to different ranges of B^, as 
shown: specifically, B^ < 0.58 in the lowest (blue) low-tan j3 region, consistent with the 
current upper limit at 95 % C.L. 
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Figure 8: The branching ratio Bgy = B{B Xg'^) x 10"^ in the [tan j3 ^Mh^) plane. The 
same CPX scenario with $^ = $3 = 90° is taken as in Fig. The different shaded 
regions correspond to different ranges of Bs-y, as shown: specifically, 3.03 < Bg^ < 4.07 in 
the upmost (blue) high-tan(5 region, consistent with the current experimentally allowed 2-a 
region, 3.03 < Bs^ < 4.07 [32]. 
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Figure 9: The ratio Rbtv in the (tan/5,M//J plane. The same CPX scenario with $^ = 
$3 = 90° is taken as in Fig. [^/or two values of p: p = 1 (left) and 10 (right). The different 
shaded regions correspond to the regions allowed at the l-a and 2-a levels by the recent 
BELLE and BABAR results: = 1.0 ±0.38 [5,35j. In the right frame, specifically, the 

2-a excluded regions are shown as Rbtv > 1-76 (in the high-tan (3 region) and Rbtv < 0.24 
(in the middle-tan j3 region). 
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Figure 10: The experimental constraints from B{Bs —>■ P'^p") (95 %), B{B — > Xg'y) (2 a), 
and Rbtv (1 <^) in the (tan/?,M//J plane for two values of p. The same CPX scenario 
with $^ = $3 = 90° is taken as in Fig. 
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Figure 11: The region allowed experimentally by the measurement B{Bd t^t~) < 4.1 x 
10~^ (90 %) [36] (upper frames) and the regions where the SUSY contribution is smaller 
than the measured values of B^-B^ mass difference [37] (middle frames) and B^-B^ mass 
difference [29] (lower frames), in the (tan/3,M/^J plane. The left three frames are for 
p = 1 and the right ones for p = 10. The same CPX scenario with $^ = $3 = 90° is taken 
as in Fig. 
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Table 4: The contents of the extended SSPARA_H(IP) 
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